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Intrinsic properties of the space itself and quantum fluctuations of its geometry are sufficient to 
provide a mechanism for the acceleration of cosmological expansion (dark energy effect). Applying 
Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy approach to self-consistent equations of one- 
loop quantum gravity, we found exact solutions that yield acceleration. The permanent creation and 
annihilation of virtual gravitons is not in exact balance because of the expansion of the Universe. 
The excess energy comes from the spontaneous process of graviton creation and is trapped by the 
background. It provides the macroscopic quantum effect of cosmic acceleration. 

PACS numbers: Dark energy 95.36.+x; Quantum gravity 04.60.-m 
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I. INTRODUCTION 



To explain the dark energy effect [H, 0a cosmological constant, hypothetical fields or some modifications of physical 
^Jy laws have been proposed (see review We show that quantum fluctuations of the metric (gravitons) and their 
back reaction on the isotropic and homogeneous (on average) background provide the mechanism for cosmological 
acceleration. The dark energy effect is a consequence of the vacuum polarization and graviton creation by non- 
| stationary gravitational field of the Universe. 

. The energy density of gravitons is a functional of the background geometry. In the non-empty Universe the 
£^ ' background geometry is defined by all contributing cosmological subsystems — by gravitons, matter and radiation. 
iy-) , The combination of conformal non-invariance with zero rest mass of gravitons (unique properties of the gravitational 
field) leads to a macroscopic quantum effect: condensation of gravitons in a quantum state with wavelength of the 
\ order of the distance to the horizon. In the process of the evolution of the Universe, the density of gravitons, as a 
result of their condensation, starts dominating over the sum total of the energy density of other subsystems of the 
| cosmological media. The self-consistent state of background and gravitons, which evolves asymptotically, represents 
self-polarized vacuum in the de Sitter space. In the subsequent paper Q we show that the regime of the de Sitter-like 
expansion is beginning to form in the current Universe which is consistent with dark energy observations 
■ i-h , We present three new exact solutions for the one-loop quantum gravity. Two of these provide cosmological expansion 
with acceleration. The de Sitter solution is one of these. All exact solutions can be found when the theory is 
presented as Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy equations for moments of the graviton 
spectral function. The same results follow from direct calculations of quantum field operator functions, state vectors 
and the graviton spectral function. 

We operate in the framework of one-loop quantum gravity because the theory cannot be renormalized in higher 
loops. Problems arising in two-loop theory are described, for example, in [6|. However, the effect of condensation of 
gravitons is created by general properties of gravitational field (conformal non-invariance and zero rest mass of gravi- 
tons). These properties will probably remain in the future comprehensive theory which includes quantum gravitation. 
The results of one-loop theory in which only the gravitational field is taken into account, are mathematically robust 
due to the finiteness of one-loop quantum gravitation [7j . In our case, the finiteness is provided by the compensation 
of diverged contributions of gravitons and ghosts to observable quantities. 
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II. EQUATIONS 

The complete classic (non-quantum) theory of back reaction of scalar, vector and tensor fluctuations on the isotropic 
space-time background has been developed in Q. In this paper, we consider the quantum theory of tensor fluctuations. 
Our model of the empty Universe consists of the background and gravitons only. In the self-consistent theory of 
gravitons, the macroscopic metric is described by regular Einstein equations 

Rl - ±&*R = Kmf* grav) +f* ghost) \V) . (1) 

The stress tensor of gravitons TH grav ^ and ghosts Th ghost ^ should be obtained by solving operator equations of motion 
and averaging over a quantum ensemble Note the average stress tensor of nontrivial ghost fields interacting with 
gravity must appear in the right hand side of ([T]) because there are no gauges that eliminate the diffeomorphism 
group degeneracy in the General Relativity. Our gauge selection was based on two principles. First, both background 
and gravitons should be considered in the same reference frame. Second, the gauge should provide automatically the 
one-loop fmiteness. We found the only gauge that satisfies both conditions comes from the set of synchronic gauges. 

Our calculations were done in the frame of one-loop approximation over quantum fields. In the flat isotropic 
Universe, the equations (fT]) read 

3H 2 = *e g = + ±Wx , -2H - '3H 2 = x Pg = -^D + , (2) 

where H = a/a is the Hubble function and a(t) is the scale factor. Here D and W\ are moments of the spectral 
distribution function of gravitons that is renormalized by ghosts. The moments are: 

D = W + 3HW a , 

k 2m ( - \ (3) 

w ™ = ^ ( EtelV&VwIp) - (gh\hh\gh)\ , m = 0, l, 2, ..,oo. 

Here and later the dots are time derivatives. Heisenberg's equations for Fourier components of the transverse 3-tensor 
graviton field and Grassman ghost field are: 



k 2 



3#6» k + — k = 



(4) 



Taking account of normalization of fields in accordance with ([2]) and ([3]), canonical commutation relations for gravitons 
and anticommutation relations for ghosts read 



a 



Km V'kv 



8x 



-ihSkk'Scrcr' , 

(5) 



Equations ([H © , Q and quantization rules §5§ have been obtained by the path integral 0, [T(| ■ They have been 
obtained from the class of synchronic gauges that automatically provide one-loop fmiteness of observables. One-loop 
effects of vacuum polarization and particle creation by background field are contained in equations ((4]) for gravitons and 
ghosts. These equations are linear in quantum fields but their coefficients depend on the non-stationary background 
metric. Correspondingly, in the background equation ([2]) we keep the average values of bilinear forms of quantum 
fields only. In this model, quantum particles interact through a common self-consistent field only. 

We would like to point out that our model contains the graviton energy stress tensor but the short wavelength 
approximation is not used. The separation of the total metric into background and gravitons is not done on the 
basis of space scales hierarchy but on the basis of symmetry criteria. The metric of the isotropic 3-space belongs to 
3-scalar representation of the group 0(3) while the gravitons are represented by transverse and traceless 3-tensor. 
The mathematically rigorous method of separating the background and the gravitons, which ensures the existence of 
the graviton energy stress tensor is based on averaging over graviton polarizations: (V'f) = if all polarizations are 
equivalent in the quantum ensemble. The properties of the theoretical model are quite specific, but it is due to this 
that we can study the effect of condensation of long wavelength gravitons in the isotropic universe. 



3 



III. BBGKY HIERARCHY AND EXACT SELF-CONSISTENT SOLUTIONS 

The following BBGKY hierarchy has been obtained from ((4]) by a standard procedure 

D + 6HD + 4Wi + I6HW1 = , 



W m +3(2m + 3)HW m + 3 (4m 2 + 12m + 6) H 2 + (2m + l)H 



w m + (6) 



+2m 



2 (2m 2 + 9m + 9) H 3 + 6(m + 2)HH + H W m + AW m+l + 8(m + 2)HW m+1 = 0, m = 1, 



In the infinite BBGKY chain ([6]), each equation connects two neighboring moments ([3]) of spectral function. The 
system of equations (J2|), ([6|) has at least three exact self-consistent solutions. Two of them are the following 

a z \ a 2/ a b J a 1 \ a 4, , , 

K 2m n K 2 n H 

W m = -24( T ir^ r ln^, m^2; H 2 =±\\n^ + ^ , 

a zm ao a z ao a 

where K 2 , ao, C are arbitrary constants. Evolution scenarios for ([7]) are simplest with C = 0. The first solution 
(upper signs in ([7])) describes the Universe that was collapsing in the infinitely remote past to the state with the 
minimal scale factor a m i n = ao, and then began to expand with acceleration a/a = K 2 /2a 2 . Asymptotically, it 
becomes logarithmically slow and reads 

a(t) ^ Ktln 1/2 (Kt/a ), t -> 00 . 

The second solution (lower signs in ((TJ) corresponds to the Universe creation from a singularity, expanding to the 
maximal scale factor a max = ao, then subsequently collapsing and ending in a final singularity. 
The third solution describes the graviton vacuum in the de Sitter space. It reads 

H= 1 -Vm, a = a e-, D = - Vr , W m+1 = - ^ + J *" + § W m , m> l. (8) 
3 2(m + 2) 

One can show that this solution is stable against small perturbations. 

From the recurrence relation for moments ([8]), we can evaluate graviton and ghost characteristic wave lengths as 
follows 

/ tj/. 1 / q 

(9) 

Quantum fluctuations of these wave lengths dominate in the formation of observables. As it is seen from ([9]), in the 
process of the Universe exponential expanding characteristic values of k rapidly shift to the region of exponentially 
large wave numbers, and all observables behave as constants with time. We will see below that this situation takes 
place if the graviton and ghost spectra are flat in k — space of conformal wave numbers, and divergent contributions 
from graviton and ghost integrals compensate each other, so values of observables are generated by finite differences 
of these integrals. 

IV. GRAVITON AND GHOST STATE VECTORS IN DE SITTER SPACE 

Below we show that the de Sitter solution follows from direct calculations of quantum field operator functions, 
state vectors and graviton spectral function, i.e. independently of BBGKY-approach. Due to the exact solution (jHJl, 
there is a vacuum state vector such that exact solutions of operator equations ([4]) become self-consistent in the de 
Sitter space averaged over this vector. These exact solutions are 




a V K (10) 
fc = ~y^[ak/(aO+0-k/*(a:)] , 
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where f(x) = (1 — i/x)e lx , x = krj, r\ = J dt/a. Substituting (fTOf into ([5]), we get commutation/anticomutaition 
relations for operator constants: 

[Ckcr, Cfc V ,l = 5kk"W , (11) 



K, Ok'] + = <$kk' , [Ac, Ac'] + = -<W ■ (12) 

In accordance with (|TT|) . the space of graviton states is constructed over the standard Fock basis |nko-)- Note that 
quantum occupation numbers nk CT = 0, 1, 2, cannot be interpreted as real graviton numbers. In non-stationary 
space quanta with arbitrary momentum are not in the mass shell. So, formally, these are virtual particles, and they 
behave as real particles asymptotically only (when kr\ 3> 1). Quanta with the characteristic wave length ([9]) are 
essentially virtual particles. Thus, occupation numbers n ka are parameters of polarized vacuum. The vacuum state 
vector is a product of superposition of states with different occupation numbers: 

ker n kcr n k(T 

where C„ ktT is the amplitude of the state with the occupation number n ka . 

According to ©, (|3]), the ghosts that correspond to the set of synchronic gauges redefine the moments of 
spectral distribution function additively. In the one-loop quantum gravity they should provide finiteness values of 
the observables. From the physical point of view, ghost fields act as compensators. We will show later that ghosts 
provide exact zero energy of quantum (and quasi-classical) gravitational waves, which cannot exist in de Sitter space. 

Grassman units are extracted from ghost operator constants multiplicatively. They read 

a k = tta k , <5k = uaj , /3 k = ubk , /3k = ub^ . (14) 

By definition, uu — —uu = 1, so operators a^, a^, 6k, 6^, which were introduced in (fT4")l . satisfy standard Bose 
commutation relations. Then, the procedure to construct the ghost space of states is obvious. The ghost state vector 
is 

1 = nz^i "k> nE^i ^> > e i^i 2 = e i^i 2 = 1 • (is) 

k «k k rik nk "k 

V. CALCULATION OF OBSERVABLES 

Further, we make use of the de Sitter solution a = —(Hrj)^ 1 and variable of integration x ~ kr\. Substituting (fTQ|) 
into §5§ and averaging over state vectors (TT3]) , (fT5|) one gets 

oo 

W ™ = ^?^ 2m+2 / ^^{^(^e)!/^ , (16) 



where 

U k (wave) = ^2(g\ci a c k a\g) - (gh\a£a k \gh) - (gh\b£b k \gh) ; 



Uk(cr) = ^E^I c ka c ik- CT l.9> - (9h\a^b\\gh) ; 
U k (ann) = ^ ^ ( g ' c -k-<xC kCT \g) - {gh\b _ k a k \gh) = U£( cr) 

a 

Here U^ wave ^ is the spectral parameter of quantum waves, which are real gravitons if krj ^ 1, and C/k( C r)i t^k(arm) 
are spectral parameters of quantum fluctuations that emerge in processes of graviton creation from the vacuum and 
graviton annihilation to the vacuum. 

Two conditions for W m , the absence of divergences and absence of t dependence, should be considered jointly. 
Because of k = x/rj dependence, U^ wave ) term in (116|) is time independent only if U^ wave ) does not depend on 
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k. If, however, U^ wave ) = const (k) ^ then this integral does not exist because of the |/(a;)| 2 — > 1 condition at 
x — > oo. This is the reason why an exact compensation of graviton and ghost contributions in U^i wave \ is a mandatory 
requirement. The compensation condition leading to U^ wave -\ = is 

I^J = M»J = |fln_J. (17) 

This result has a simple physical interpretation. Quantum gravitational waves, whose equation of state is not p = —e, 
cannot exist in de Sitter space with the self-consistent geometry. 

Analogously, spectral parameters t/k(cr) = U, Ukiann) = U* cannot depend on k also. But corresponding integrals 
in l|16p do not lead to obligatory divergences because oscillating functions ~ e ± 2M are integrated at x — ► oo. Such 
integrals may be redefined to have certain finite values. Thus, in (|16p we have flat graviton and ghost spectra, 
Uk( wave ) = 0, f7 k(cr) = U* (ann) = U and 




(18) 



I <An | — | B n | — y*Pn 



where V n is any universal distribution and n = nVin) is the average number of gravitons with wave lengths that 

n=0 

are near the characteristic value ([5]). 

From (fT7|) and (fT8| it follows that U = if all amplitudes are real. We have U ^ if some amplitude are complex 
and phase dependencies of occupation numbers are different in graviton and ghost sectors. Observable quantities are 
proportional to U + U* = 2N g , where N g is a generalized parameter of graviton vacuum state. To use some model 
distributions (Poisson distribution, for example) one can see that N g ~ n. 

Integrals in (fT6)l should be redefined. They can be calculated in the following way 



lim [ dxx 2m±1 e-^- 2l)x = =F(-1)' 



(2m ± 1) 



( ^ J 2 2m + 1±1 ' 

° oo (19) 
2i lim [ dxx 2m e-«- 2 ^ = (-l) m+1 ^1 . 



<-> o J v ' 2 2r - 



The physical interpretation of (119|) is as follows. In any instant of time, the procedure of redefining the integrals in 
(fT9)) selects the contribution of virtual gravitons with wave lengths that are near characteristic value ([9]) only, and 
eliminates contributions of all other gravitons. This redefining procedure provides the existence of the exact solution 
© of BBGKY hierarchy. 

Zero moment Wq is not included in the equations of the theory. Therefore, to avoid a logarithmic infrared singularity, 
the function under the integral sign in (fT6|) is differentiated initially, then derivatives are combined as D = Wq+SHWq, 
and redefined integrals are calculated. The results are 

D = -^^H\ w m ^t^f^ { 2 m ~l)\(2 m +l)( m + 2)x^^H 2m+2 , m > 1. (20) 

It is easy to verify that expressions (|20| satisfy to recurrence relations (JSJ, following from the BBGKY hierarchy 
immediately. 

In accordance with ^ and (|2H)) . the energy density and the pressure of virtual gravitons (that generate the back 
reaction of quantum fluctuations on background geometry) are 



One can also obtain the equation (|2ip by other methods. The vacuum amplitude of graviton field in the one-loop 
approximation can be presented via a path integral, which, as it is known [111 ], can be calculated exactly in the de 
Sitter space. One can calculate this path integral under boundary conditions that correspond to the state vector (fT3j) . 
(p~5|> . under the assumption that the de Sitter geometry is self-consistent. A variation of result of this calculation 
over the scale factor leads to the convolution of Einstein equations |T]) and thus to the expression (e g — 3p g )/4 = e g 
that coincides with (|2"Tj) . The same result can be obtained by Schwinger-DeWitt method using the effective action 
formalism. 
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VI. ENERGY BALANCE 



Equation (|21[) is a consequence of the balance equation for the energy that is emerging to the space due to graviton 
creation and vanishing due to graviton annihilation. The characteristic energy of gravitons in these processes is 
hw r*> HH. Total probabilities of graviton creation and annihilation processes (normalized to the unit volume) w cr 
and w ann are proportional to the phase volume of one graviton lu 3 /Sir 2 ~ H 3 /3ir 2 . The exponent of the background- 
graviton coupling constant is unity if lo ~ H . Thus, we obtain for w cr and w ann the following estimations 

w °r = A ff3 (^k + l)(iV-k+l), w ann = -^H 3 iV k AU k . (22) 

where a = 0(1), N±k ~ N g /2 is the average number of gravitons with wave lengths that are near characteristic value 
©. Finally, we get the balance equation in the form 

e g = hw (w cr - w ann ) = "hH*(N k + AU k + 1) ~ ^-hN g H A . (23) 

This estimate with accuracy of numerical factor of the order of unity coincides with (|2~lj) . From (|2"2"|) and (j2"3"| one can 
see that the non-zero effect occurs because of quantum spontaneous process of the particle creation. 

If non-gravitational contributions to the vacuum energy are compensated exactly (due to, e.g., a consequence of 
Supersymmetry) , it follows from ^ and (|2"Tj) 

IP = . (24) 

According (|24p. iJ 2 ~ ft" 1 , and it clearly demonstrates that the acceleration of the Universe expansion is a macroscopic 
quantum effect in the graviton vacuum at the scale of the Universe (analogous to such known macroscopic quantum 
effects as super- fluidity and super-conductivity). Substituting the observable value of Hubble constant into (|24[) . we 
get an estimate of the total number of virtual gravitons under the de Sitter horizon N g ~ 10 123 . This could be a 
possible interpretation of the s DE /epi an k ~ 1CV 123 ratio. 



VII. CONCLUSION 



Virtual gravitons with wavelength of the order of the horizon must appear and disappear in the graviton vacuum 
because of massless and conformal non-invariance of the graviton field. A non-zero balance of energy is due to the 
pure quantum process of spontaneous graviton creation, in other words, due to the uncertainty relation. Numerical 
integration of the system of equations (O shows several types of solutions with acceleration. In this paper, 
we considered the physics of this effect by the example of the exact de Sitter solution, which is the only one that 
corresponds to the state of equilibrium of graviton vacuum, and this is the reason why this exact solution exists. 
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